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Abstract 

We introduce an unusual Mexican hat potential, a piecewise parabolic one, 
and we show that its vortex solutions can be found analytically, in contrast 
to the case of the standard \^\ 4 field theory. 



Spontaneous symmetry breaking has been studied traditionally through the 
\1/| 4 potential, 



This potential is usually called a Mexican hat potential by particle physicists, 
even though the actual Mexican hat has somewhat different wings. The above 
potential is relevant not only to particle physics, but also to condensed matter 
physics. In superffuid helium II, for example, it is used for writing down the free 
energy, and it leads to a lot of interesting physics. 

Of particular interest are the vortex solutions that this potential admits. In 
the context of superfluid helium, these were studied quite a long time ago[l]. Due 
to their inherent nonlinearity, the corresponding field equations have to be solved 
numerically. There are analytic approximations for the vortex solutions[2], but 
the exact solutions can only be found numerically [3]. 

Vortices play an important role in many fields of physics. It may be of interest 
then to have analytic vortex solutions, that could be used in further calculations 
involving vortices. Such analytic solutions would also be of pedagogical inter- 
est, since they would show explicitly the various properties of a generic vortex 
solution. 

In this paper we study an unusual Mexican hat potential that admits vortex 
solutions. These vortex solutions are analytic though, and are expressed in terms 
of standard special functions of mathematical physics. Not only can one study 
therefore the properties of the solutions analytically, but one could also use them 
in other calculations involving vortices. 

Our potential is 



This potential is a Mexican hat potential, but with a kink at ^ = 0. It is a piece- 
wise parabolic potential. Such potentials, but with real order parameters, have 
also been used in Ginzburg-Landau theories of oil-water-surfactant mixtures[4], 



V(V) = A(|^| 2 -t; 2 ) 2 . 
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because they lead to analytically solvable equations. 



The static free energy is then 
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If we measure \& in units of v, x, y and z in units of I/7, and F in units of Xv 2 /^/ 3 , 
where 7 2 = 2m\/h 2 , then we obtain the dimensionless free energy 
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This is minimised when 
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This is the basic field equation. We note first that it admits one-dimensional 



solitonic solutions: 



tf(s) = .i_fo (i_ e -l*-* ol 
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zq being arbitrary. Both * and <9\I/ / dz are continuous at z = zq in this solution. 
It is worthwhile to note that the one-dimensional solitonic solutions of the I*) 4 
theory are also known analytically, unlike the |\l/| 4 vortex solutions. 

Let us now concentrate on the solutions of our potential of Eq. 2 for a single 
isolated vortex. We are thus looking for solutions 



(7) 



where x = rcos0, y = rsin0, n being a positive integer, with ip(r) > 0. In that 
case Eq. 5 reduces to 



d 2 ip 1 d^/j 



,n 



+ = -1. 



(8) 



dr 2 r dr r 

Note that this equation is a linear differential equation, albeit an inhomogeneous 
one, while the corresponding one in the |\I/| 4 theory is nonlinear. The dimen- 
sionless free energy for a cylinder of radius R, and length L along the z-axis, 
becomes 

' n " "> ■ ' ' — 1 (9) 
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Eq. 8 clearly shows that ip — > 1 as r — > oo. In particular, ip — (n 2 /r 2 ) in that 
limit. Hence the free energy contains a piece proportional to n 2 that diverges 
logarithmically as R — > oo. For large R then it is clear that the solution n — 1 is 
best, as is usually the case for isolated vortices. 

It is interesting to examine the behaviour of ijj for small r. Eq. 8 easily yields 
that ip is linear in r for small r, if n — 1, while ^ mr 2 / (n 2 — 4) if n > 2. We have 
then an important difference between our potential and the |\l/| 4 theory when it 
comes to highly quantized vortices. Namely, while ip oc r n in the |^| 4 theory, 
here ip oc r 2 for n > 2. 

Let us now multiply Eq. 8 by rip and then integrate from to R. We shall 
find a relation that enables us to write the free energy in the exact form 
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This equation is valid for all n. 

We now present the solutions of Eq. 8, beginning with the case where n is a 
positive odd integer. Let us define first the function 

g(r) = l 7 / e- rcose cosn^^, (11) 

yv 1 sin(rar/2) Jo v 1 

with n = 1,3,5... Note that g — > 1 as r — > oo, and that g(0) = 0. In particular, 
when r is small, g ~ 7rr/4 for n = 1, and g ps r 2 /(n 2 — 4) for n = 3,5,7... We can 

easily prove the identity 

n 2 n f w / 2 d 



r —\e- rcose (sm6cosn6 + -smn6)]d6. (12) 
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Its r.h.s. is equal to 
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Eqs. 12 and 13 imply then that g satisfies Eq. 8. It is thus a particular solu- 
tion of Eq. 8. The most general solution is g(r)+CiK n (r)+C2l n (r), in terms of 



modified Bessel functions. However, K n diverges at the origin and I n diverges at 
infinity,while g is well behaved at both limits. Thus we must have c\ = c<i = 0. 
Hence 

*p(r)=g(r). (14) 

In particular, for n — 1 we get 

f7 r/2 
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e- rcos9 cos6d6]. (15) 
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This is the exact solution representing a singly quantized vortex, and the main 
result of this paper. It can be written in the form 

* = |e i *[/ 1 (r)-L 1 (r)], (16) 

where ii(r) is a modified Bessel function and L^r) a modified Struve function[5]. 
The corresponding ip{r) is shown in Fig. 1. 

Actually, one expects this result, because the function 2ifj(—ir)/-K can be 
shown to satisfy the standard Struve differential equation [5]. It is thus, in fact, 
that one gets the idea of using the function of Eq. 11, since the integral rep- 
resentations of the Struve functions involve integrals such as those used in Eq. 
11. 

We can easily show, using Eq. 15, that ip(r) « irr/4 for r — > 0, for this singly 
quantized vortex. Furthermore, the asymptotic properties of the modified Struve 
function give, for r — > oo, 

Thus the behaviour of our singly quantized vortex at infinity and at the origin 
resembles that of the corresponding |\1/| 4 vortex in these limits. Indeed, the ip for 
the |\1/| 4 vortex is linear at the origin and tends to 1 — (2r 2 ) -1 at infinity. 

Let us now calculate the free energy, using eq. 10. The surface term is 2/i? 2 , 
for large R. We can use furthermore the identity [5] 

r[l - 7r(/i - Lx)/2] = (L - / )] + \{h ~ U) (18) 



to get 

^^^(WX + l/Vo-LoM, (19) 

But we also have[5] the properties L (0) = 0, J o (0) = 1, L (R)-I (R) &-2(R~ 1 + 
R~ 3 )/n, and /<f (I -L )dr « [2 ln2 J R+2 7 - J R- 2 ]/7r, where 7 =0.5772157 is Euler's 
constant. Using all these, we obtain the final result for the singly quantized 
isolated vortex: 

ln#+( 7 + ln2-l) + -]-^ln(i.3LR). (20) 



2nL w ' 2R 2 

The corresponding result for the |\l/| 4 singly quantized vortex[l], for a cylinder 
with dimensionless radius R, is ln(1.46i?). 

For vortices with circulation n = 3,5,7,..., we can find the free energy by 
combining eqs. 10 and 14, in which case we get 

~ n 
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In order to evaluate this integral, we split it into two pieces, one piece from 
to uj, and a second piece from u to it/2, where a; is a small number such that 
sincj « uj and e~ Rsm0J rs 0. We can then drop the exponential in the fj^ 2 piece, 
since R is large, and we can let sin 9 ~ 9 in the Jq piece. Thus 



r5l- t -»-*-«] + »f (22) 
Jo 9 Juj sin^ 
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Now the integrals can be done easily, leading to 

, (n-l)/2 2 1 

^«„'(b,fl +7+ h.2- g sn--) (23) 

for n = 3,5,7.... For example, f/2nL w 9 ln(0.345i?) for n = 3, while the corre- 
sponding |\I/| 4 vortexfl] has a free energy f /2-kL rj 91n(0.38i?) for a cylinder of 
length L and dimensionless radius R. For n = 5, 7, 9 the energies of our vortices 
are 25^(0.20^), 49^(0.1447?) and 81 ln(0.11.R). Of course, the case n = 1 is 
the physically interesting one, since it is energetically preferred. We shall find 
nonetheless the vortex solutions for the case when n is even as well, since the 



results are simple and analytic. 

We examine now then the case n=even, which is actually simpler than the 
case of odd n. We can easily show that the function 

h ( r ) = n £ (^-^- 1 ) ! r 2fc-n ( _ 1)fc+ f 2 n-2 fc -l ^ 

is a particular solution of Eq. 8, if n is even. The most general solution would be 
•0(r) = h(r)+ciK n (r)+c 2 I n (r). We must have c 2 =0, otherwise the solution will 
diverge at infinity. We also note that the most divergent part of K n (r) at the 
origin is (n — l)\2 n ~ 1 r~ n , while the most divergent part of h(r) at the origin is 
n\r~ n (— i) n / 2 2 n_1 . The solution must not have a r~ n piece at the origin, hence we 
must have C\ = — n(— l) n / 2 . Remarkably, this choice ensures that all divergences 
at the origin disappear, leaving us with i/)(0) = 0, as expected for a vortex. Thus 
the exact vortex solution for n=even is 

1>(r)=h(r)-n(-l) n ' 2 K n (r). (25) 

For example, the n = 2 vortex has 

*P(r) = l-^ + 2K 2 (r), (26) 

while the n = 4 vortex has 

^(r) = l-^ + -£--4K- 4 (r). (27) 

We can now evaluate / /27rL using Eq. 10. For the vortex with n — 2, straight- 
forward integration using the solution of Eq. 26 yields, for large R, 

f 8 32 

^ = [41nr + 2tK 1 {t) + 4K (r)}« + ^ - ^ ^ 

which reduces to 

' =41 nj R + 4 7 -41n2-2 + J- - ^ ^ 4 ln(0.54i?). (29) 
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The corresponding result [1] for the |\1/| 4 vortex with n = 2 is 41n(0.59i?). We 
note that here too, just as with the cases n — 1 and n — 3, the coefficients of 



R within the logarithm are quite close to those that are appropriate for the |^| 4 
vortices. In detail, our vortices have the coefficients 1.31, 0.54, 0.345 for n — 1, 
2, 3 respectively, while the |\I/| 4 vortices have 1.46, 0.59, 0.38. 

Let us calculate the free energy of the other n=even vortices as well, using 
Eqs. 10 and 25. The integral J R rK n (r)dr is needed. We can easily show that 

—— i 

±[-rK n _ 1 (r) + 2nJ2(-iyK n - 2t (r)+n(-l) n / 2 K (r)} = rK n (r), (30) 
ar i=i 

whence 

— — i 

_Z_^_ n + n 2[ lnjR + 7 _ ln2 _ (31) 

The divergences that the various modified Bessel functions of Eq. 30 exhibit at 
the origin cancel, confirming thus the correct evaluation of the free energy. 

We have thus completed the evaluation of the exact isolated vortex solutions 
for our unusual Mexican hat potential. We have been able to obtain analytic 
solutions for any circulation n, odd or even. This success is due to the fact that 
the field equation is not nonlinear. It is a linear inhomogeneous equation, and it 
can therefore be solved exactly. We find the usual features that a vortex must 
have, i.e. that ip vanishes at the origin, while it tends to a constant at infinity. 
The free energy depends logarithmically on the radius of the vortex, as expected. 
Finally, the actual values of the free energy are quite close to those of the free 
energy for the standard l^j 4 vortex. 

Aside from the pedagogical value of these solutions, especially when discussing 
spontaneous symmetry breaking, our results may be used in more elaborate cal- 
culations involving vortices, in many fields of physics. 
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Figure Captions 
Figure 1: The solution for a singly quantized isolated vortex. 



